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Abstract 

Constitutive equations are derived for enthalpy recovery in polymeric glasses after 
thermal jumps. The model is based on the theory of cooperative relaxation in a version 
of the trapping concept. It is demonstrated that some critical temperature Tcr and some 
critical degree of crystallinity fcr exist in a semicrystalline polymer above which the 
energy landscape becomes homogeneous and structural relaxation ceases. 

1 Introduction 

This note is concerned with the kinetics of enthalpy recovery in polymeric glasses after 
thermal jumps. Physical aging of polymers has been the focus of attention in the past three 
decades, see monographs 0, ||, 0, ^, ^ and review articles [§, |^, ||. Despite significant 
successes in the study of out-of-equilibrium dynamics in disordered materials achieved in 
the past decade (see reviews ||^, |10|, |TT|, [13], ^^), it is difficult to mention a model that 
adequately predicts structural relaxation in polymers. Even the mode-coupling theory |]15 



(the most advanced among molecular models) fails to describe slowing down in the response 



of amorphous polymers below the glass transition temperature Tg [0, [T6|, |T3- As a reason for 



this conclusion, the neglect of cooperativity in the molecular reorientation is conventionally 



mentioned |Tq, 19 



Enthalpy relaxation provides one of the most convenient ways to study aging of glassy 
polymers [0, |2y, 0, ^ |23|. The kinetics of enthalpy recovery is in close connections with 
time- varying changes in static and dynamic moduli ^|, specific volume ^ 



and yield stress p9|, |3T|, as well as with transition from ductile to brittle fracture 
31]. For some polymers (e.g., polyesters), DSC measurements demonstrate a dramatic 



effect of the waiting time whereas the free volume fraction determined by PALS shows 



no pronounced decrease induced by physical aging p6| . 

Constitutive models for structural relaxation in glassy polymers are conventionally con- 
fined to amorphous media (as exceptions, we would mention Refs. |]37| , ^). This may be 



explained by a belief that crystalline lamellas do not affect the process of structural relax- 
ation occurring in amorphous regions This hypothesis may be accepted provided that 



only one kind of amorphous domains exists in a semicrystalline polymer. A number of ob- 
servations demonstrate, however, a more complicated picture, where two different kinds of 
amorphous regions co-exist in a partially crystallized polymer: amorphous layers between 
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lamellas and amorphous domains between spherulites This results in a 



question whether the rate of physical aging is a linear function of the degree of crystallinity 
which vanishes when the degree of crystallization / reaches 100% only (as it is postulated in 
Fig. 6 of Ref. [^), or the rates of structural relaxation substantially differ in interlamellar 
and intralamellar zones which implies that physical aging may decline even in a partially 
crystallized material. 

The objective of this study is to derive a constitutive model for enthalpy recovery which 
is based on the theory of cooperative relaxation in a version of the trapping concept |TB|, H 



45| , To simplify the analysis, we do not distinguish explicitly between interlamellar 

and intralamellar regions, but presume that the "averaged" distribution of energies of traps 
may depend on the degree of crystallinity. An amorphous region is treated as an ensemble 
of mutually independent cooperatively rearranged regions (CRR). A CRR is thought of as a 
globule consisting of scores of strands of long chains . The characteristic length of a CRR 
in the vicinity of Tg amounts to several nanometers [Q. In the phase space, a CRR is treated 
as a point trapped in its potential well (cage). At random times, the unit hops to higher 
energy levels as it is thermally agitated. When the energy of thermal fluctuation exceeds the 
height of a barrier between cages, a CRR can change its trap. Introducing several hypotheses 
regarding the kinetics of hops from one potential well to another, we develop a nonlinear 
parabolic equation for the probability density of traps with various potential energies. This 
equation is applied to the analysis of enthalpy recovery in amorphous and semicrystalline 
poly (ethylene terephthalate) (PET). The following conclusions are drawn: 

1. the degree of crystallinity of a semicrystalline polymer affects the model's param- 
eters in a similar way as temperature influences those for an amorphous polymer (a 
crystallinity-temperature principle analogous to the time-temperature principle in lin- 
ear viscoelasticity p9|); 



2. a critical level of crystallinity < 100% exists at which the energy landscape becomes 
homogeneous. This level may be treated as an analog to the critical temperature 
for structural relaxation in amorphous polymers |^ predicted by the mode-coupling 
theory. 

The exposition is organized as follows. Constitutive equations for the kinetics of enthalpy 
relaxation are developed in Section 2. These relations are verified in Section 3 by comparison 
with experimental data. Some concluding remarks are formulated in Section 4. 



2 A model for enthalpy relaxation 



Denote by w the energy of a potential well with respect to some reference state [^]. It is 
assumed that w > for any trap and w = for the reference state. At random times, a 
CRR hops to higher energy levels in its potential well as it is thermally activated. Denote 
by q{(jj)duj the probability to reach (in a hop) the energy level that exceeds the bottom level 



of its potential well by a value located in the interval [uj^uj + duj]. Referring to [T^, we set 



q{uj) = Aexp(— Atu), where A is a material constant. The probability for a CRR in a trap 
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with potential energy w to reach the reference state in an arbitrary hop is given by 



poo 

Q{w) = / q{uj)<ko = exp(— y4w). 

J w 



The average rate of hops in a cage 7 is determined by the current temperature T only, 
7 = 7(T). The rate of rearrangement P equals the product of the rate of hops 7 by the 
probability Q to reach the reference state in a hop, 



P{w) = •y exp{—Aw). 



Denote by X the (time-uniform) concentration of traps per unit mass, and by p{t, w) the 
current probability density of traps with potential energy w. The number of relaxing regions 
(per unit mass) trapped in cages with the energy belonging to the interval [w,w + dw] and 
rearranged during the interval of time [t, t + dt] is XP{w)p{t, w)dwdt. Unlike Refs. 
we assume that not all flow units change their traps when they reach the reference energy 
level and denote by F{t,w) the ratio of the number of relaxing regions returning to their 
traps to the number of those reaching the reference state. The number of relaxing regions 
leaving their cages (with the energy located within the interval [w,w + dw]) per unit mass 
and unit time is given by X{1 — F)Ppdw. The exchange of flow units is assumed to occur 
only between the nearest neighbors on the energy landscape, that is between a trap with 
the energy [w, w + dw] and traps with the energies [w — dw, w] and [w + dw, w + 2dw]. The 
balance law for the number of flow units trapped in cages with the energy belonging to the 
interval [w,w + dw] reads 



dp 
dt 



= -(1 - F)Pp + 2 [(1 - P)Pp 



+ 2 



(1 - F)Pp 



where the subscript indices "— " and "+" refer to appropriate quantities for the intervals 
dw,w] and \w + dw,w + 2dw]. Expanding the right-hand side of this equality into 



\w 



the Taylor series, using Eq. (1) and introducing the notation T = ^'jdw'^, we arrive at the 
differential equation for diffusion over the energy landscape 



dp 
di 



52 



dw'^ 



[1 — F) exp{—Aw)p 



We adopt the Metropolis transition rates |54], 



F{t,w) 



1, P{t,w) < Pooiw), 

exp[-e(p - Poo)], p(i, w) > Poc{w), 



(3) 



An 

m is that 

under condition (3) it does not impose restrictions on the equilibrium density of trap Poc,{w). 
Referring to |ff6l, we suppose that the inequality 



where Poo{w) is the equilibrium density of traps and e > is a material parameter, 
important advantage of Eq. (2) compared to relationships suggested in Refs. 



.0 

/ p{t, w)dw -C 1 
J —00 



(4) 
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is satisfied for any t > and describe tlie initial distribution, po{w) 
distribution, poo{w), by tlie Gaussian formulas 



and the equilibrium 



Po{w} 



/27rEo 



exp 



[w 



wy 



2Sg 



1 



Poo[UJ} 



exp 



[w 



(5) 



where W, Sq, Sqo are adjustable parameters. Equations (5) imply that the average equi- 
librium energies of traps are temperature-independent (the same value W is employed for 
the initial and equilibrium distribution functions), whereas their variances strongly depend 
on T. The first assertion is fairly well confirmed by experimental data in mechanical tests 



57, 58 1, whereas the other hypothesis is in agreement with the conventional scenario for the 



growth in the ruggedness of the energy landscape with a decrease in temperature |jTo|, ^ . 
The level of disorder in an ensemble of CRRs is characterized by the configurational 



entropy per rearranging region ^0 



POO 

s{t) = —kB I p{t,w)\np{t,w)dw^ 
Jo 



where ks is Boltzmann's constant. The configurational enthalpy per unit cage h{t) is ex- 
pressed in terms of the configurational entropy s{t) by means of the conventional formula 

^ = T 
ds 

This equality is integrated for a standard one-step thermal test. 



T(t)=To (T<0), T{t)=T (T>0), 

and an explicit formula is found for the enthalpy per unit mass H 
enthalpy per unit mass AH(t) = H(t) — H{0) is given by 



(6) 

Xh. The relaxing 



AH{t) = A / po{w)\npo{w) — p{t,w)\np{t,w] 
Jo L 



dw (7) 



with A = ksTX. Introducing the dimensionless variables w = Aw and t = t/to, where to 
is the characteristic time of aging, and setting T = A^Tto, W = AW and Sfc = AE^, we 
arrive at the constitutive model, Eqs. (2), (3), (5) and (7), with six adjustable parameters 
W, Sq; Soo, r, e and A. These kinetic equations substantially differ from conventional 
relations for enthalpy recovery in polymers [^, ^ ^ because they do not refer to 
a (purely phenomenological) concept of internal clock. As an analog of the material time 
r in Eqs. (2), (3), (5) and (7), we may mention the parameter e that characterizes the 
rate of changes in the energy landscape. An important difference between e and r is that 
e is independent of the current energy landscape (but, in general, temperature-dependent), 
whereas the parameter r is conventionally expressed in terms of the current enthalpy H by 
means of the Narayanaswamy or the Adam-Gibbs equations. 



3 Comparison with experiments 

It is easy to check that the quantities W and F are interrelated: when one of them is 
fixed, the other may be chosen to characterize the time scale. For convenience of numerical 
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simulation, we fix W and determine F by fitting observations. The value W = 2.5 ensures 
that inequality (4) is satisfied with a high level of accuracy. 



We begin with experimental data for two amorphous PETs exposed in Ref. [^]. First, 
observations are fitted for a homopolymer to determine the parameters Eq? ^oo, T, e and A 
that ensure the best approximation of observations. Afterwards, the amounts Sq, Sqo and 
F are fixed, and the quantities e and A are determined by matching experimental data for 
a copolymer. Figure lA demonstrates fair agreement between observations and results of 
numerical simulation. Using the specific gravity g = 1.34 g/cm^ at T = 65 °C we find 
the volume concentration of relaxing regions in PET, E = 6.39 ■ 10^6 m-^ which is rather 
close to S = 4.1-10^^ m~^ for polycarbonate and H = 7.7-10^^ m~^ for poly (vinyl acetate) 
Il65| found by PALS. To ensure that adjustable parameters are rather robust with respect to 
changes in material properties and conditions of the test (a decrease in Tg by 5 K and an 
increase in Tq by 15 K), we repeat this procedure using experimental data presented in Ref. 
||66|| . Figure IB demonstrates fair agreement between observations and results of numerical 
simulation for similar values of Eg, e and A. An increase in F depicted in Figure IB compared 
to Figure lA seems quite natural, because this parameter is determined by the difference 
between the glass transition temperature and the annealing temperature, AT = Tg — T. A 
decrease in Sqo may be associated with a drop in Tg, because it is conventionally accepted 
that the glass transition temperature grows when molecular mobility becomes more restricted 



(e.g., because of an increase in the number of crosslinks [|30[ or in the molecular weight 
which may be associated with the growth of inhomogeneity of the energy landscape. 

To study the effect of the annealing temperature T on the kinetics of structural recovery. 



we approximate experimental data for amorphous PET obtained in Ref. Figure 2A 

demonstrates fair agreement between observations and results of numerical simulation in 
the temperature range from Tg — 25 to Tg — 5 °C. The standard deviation of energies of 
traps in thermodynamic equilibrium is plotted versus the degree of undercooling AT in 
Figure 2B. This figure demonstrate that the dependence Soo(T) is fairly well approximated 
by the linear function 

Soo = ao + ciiAT, (8) 

where are adjustable parameters. In follows from Eq. (8) that some critical temperature 
Tcr exists at which Eoo vanishes and the energy landscape becomes homogeneous. Results of 
simulation imply that T^r = Tg + 34.11 °C which is in accord with the values of T^ found by 
fitting date in mechanical tests [Q. The quantities F and e are depicted in Figure 3A. This 
figure show that the dependences f (T) and e(T) are correctly approximated by the "linear" 
functions 

log f = 6o - &i AT, log e = Co + ci AT (9) 

with adjustable parameters bk and Ck- In the vicinity of the glass transition temperature, 
the apparent activation energy AE is calculated as 

Ai? = -i?-^^"^ 



^d{l/T) 



t=tJ 



where R is the gas constant. This equality together with Eq. (9) results in the formula 
AE = RT%i. According to Fi gure 3A, AE = 156.1 kJ/mol, which is quite comparable 



with AE = 206 kJ/mol determined for PET in mechanical tests An increase in the 
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parameter f with temperature T is in agreement with the theory of thermally activated 
processes. A surprising result is a decrease in e with temperature. It follows from Eq. (3) 
that this decrease provides an additional source for slowing down of the aging process in the 
close vicinity of Tg which has not been accounted for in previous studies. 

To assess the effect of crystallinity on the rate of structural relaxation, we approximate 
experimental data for enthalpy recovery in semicrystalline PET with various degrees of 
crystallinity /. For a detailed description of specimens and the experimental procedure, 



see Ref. |^3|]. First, we approximate observations for a sample with the smallest value of 
/ and determine adjustable parameters of the model using the steepest-descent procedure. 
Afterwards, we fix the quantities So and A and repeat matching experimental data with 
three adjustable parameters, Soo, f and e. Figure 4A demonstrates fair agreement between 
observations and predictions of the model. The equilibrium standard deviation of energies 
of traps Soo is depicted versus the degree of crystallinity / in Figure 4B. The dependence 
Soo(/) is correctly approximated by the linear function 

Soo = ao-«i/ (10) 

with adjustable parameters a^. Comparing Figures 2B and 4B, we conclude that an increase 
in / for a semicrystalline polymer affects the equilibrium distribution of cages in a way 
similar to that in which the annealing temperature T influences the energy landscape for an 
amorphous material. By analogy with an amorphous polymer, one may define the critical 
degree of crystallinity for a semicrystalline medium as the percentage of crystallites at 
which the parameter Sqo vanishes and the energy landscape becomes homogeneous. In 



contrast to the conclusions of Ref. [^, Figure 4B demonstrates that is essentially 
less than 100%. The parameters T and e are plotted versus the degree of crystallinity in 
Figure 3B. This figure shows that the dependences r(/) and e(/) are correctly approximated 
by the "linear" functions 

logr = 6o + &i/, loge = co-ci/ (11) 

with adjustable parameters and c^. 

Two kinds of crystallization of polymers are conventionally studied: thermal crystalliza- 
tion, when a specimen is annealed at a fixed temperature in the rubbery region for a given 
time (this procedure was employed in Ref. [0) and mechanically induced crystallization. 



when a sample is stretched slightly above Tg to a given extension ratio A. To demonstrate 
that the effect of the degree of crystallinity on the kinetics of physical aging weakly depends 
on the crystallization procedure, we approximate experimental data for semicrystalline PET 
samples where the level of crystallization is established by uniaxial stretching at an elevated 
temperature. A detailed description of specimens and the experimental procedure can be 



found in Ref. |^ . The extension ratio A is transformed into the level of crystallinity using 



Figure 1 of Ref. |70]. Unlike previous sets of experimental data, where e is treated as a 
function of the level of crystallinity /, we fit observations at various elongations A with the 
same value of e found for an undeformed specimen. Figure 5A demonstrates good agreement 
between experimental data and predictions of the model. The functions Soo(/) and r(/) are 
plotted in Figure 5B. This figure demonstrates that Eqs. (10) and (11) adequately describe 
observations. The parameter bi (that characterizes the influence of crystallinity on the rate 
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of rearrangement) acquires similar values for thermally and mechanically crystallized sam- 
ples. The same assertion is true for the critical degree of crystallinity fcr which adopts the 
value 33.9% for the stretched PET versus 42.7% for the PET annealed above Tg. On the 
contrary, the parameter for the PET obtained by thermal crystallization exceeds that 
for the hot-drawn PET by twice. This may be explain by the fact that stretching a rubbery 
polymer establishes an additional order caused by partial destruction of a polymeric network 
and alignment of long chains along the axis of loading |]32| . 

4 Concludung remarks 

Constitutive equations have been derived for enthalpy relaxation in glassy polymers after 
thermal jumps. The model is based on the trapping concept which treats a disordered 
medium as an ensemble of flow units rearranged at random times as they are thermally 
activated. Adjustable parameters are found by fitting experimental data for amorphous and 
semicrystalline poly (ethylene terephthalate). The following conclusions are drawn: 

1. the constitutive equations correctly describe the kinetics of structural relaxation in the 
sub-Tg region. 

2. the model predicts the existence of some critical temperature at which the energy 
landscape becomes homogeneous. 

3. for semicrystalline polymers, the model implies the existence of some critical degree of 
crystallinity at which the energy landscape becomes homogeneous. The inequality 
/cr < 100% may serve as an indirect confirmation that the processes of structural 
relaxation in interlamellar and intralamellar regions differ from one another. 

4. a correspondence may be established between the effects of temperature and the level 
of crystallinity on the kinetics of enthalpy recovery. 
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Figure 1: A — The relaxation enthalpy /\H J/g versus time t days for amorphous PET 
(Tg = 80 °C) quenched from Tq = 80 to T = 65 °C. Circles: experimental data [ P3|] . 
Solid lines: prediction of the model with Sq = 0.19, Sqo = 0.95 and f = 8.0 days~^. 
Curve 1: homopolymer, e = 0.2, A = 4.0 J/g; curve 2: copolymer (3.5 mole % of cyclohexane 
dimethanol), e = 0.17, A = 4.4 J/g. B — The relaxation enthalpy AH J/g versus time t h for 
amorphous PET (Tg = 75 °C) quenched from Tq = 95 to T = 65 °C. Circles: experimental 
data [0]. Solid line: prediction of the model with Sq = 0.15, Sqo = 0.38, f = 1.24 h~\ 
e = 0.2 and A = 4.2 J/g 





30.0 



Figure 2: A — The relaxation enthalpy AH J/g versus time t h for amorphous PET (Tg = 
76 °C) quenched from Tq = 277 to the test temperature T °C. Symbols: experimental data 
||35| . Sohd lines: prediction of the model with Eq = 0.25 and A(Tg) = 2.9 J/g. Curve 1: 
T = 71; curve 2: T = 67; curve 3: T = 60; curve 4: T = 52. B — The parameter versus 
AT °C. Circles: treatment of observations. Solid line: approximation of the experimental 
data by Eq. (8) with = 0.433 and ai = 0.013 
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Figure 3: The parameters T h^^ (curve 1) and e (curve 2) versus AT °C for amorphous PET 
(A) and versus f% for semicrystaUine PET (B). Symbols: treatment of observations p5| , 
Sohd hues: approximation of the experimental data by Eqs. (9) and (11). A: Bq = 3.686, 
6i = 0.154, Co = -2.424, ci = 0.073; B: bo = 0.618, h = 0.038, cq = 0.247, ci = 0.071 
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Figure 4: A — The relaxation enthalpy AH J/g versus time t h for semicrystaUine PET 
(Tg = 80 °C) quenched from Tq = 100 to T = 65 °C. Symbols: experimental data |Q. Solid 



lines: prediction of the model with Sq = 0-2 and A = 1.5. Curve 1: / = 14%; curve 2: 
/ = 21%; curve 3: / = 27%; curve 4: / = 32%. B — The parameter Soo versus /%. Circles: 
treatment of observations. Solid line: approximation of the experimental data by Eq. (10) 
with ao = 1.519 and ai = 0.036 
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Figure 5: A — The relaxation enthalpy AH J/g versus time t days for semicrystalline PET 
(Tg = 77 °C) uniaxially stretched at To = 100 °C and quenched to T = 40 °C. Symbols: 
experimental data [^. Solid lines: prediction of the model with Eq = 0.12, e = 0.2 and 
A = 3.38. Curve 1: A = LO; curve 2: A = 2.0; curve 3: A = 3.0; curve 4: A = 4.0. B— The 
dimensionless parameter (curve 1) and the parameter T days^^ (curve 2) versus /%. 
Circles: treatment of observations. Solid lines: approximation of the experimental data by 
Eqs. (10) and (11) with = 0.603, ai = 0.018 and bo = -0.868, bi = 0.027 
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